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Empirical mode decomposition is an adaptive signal processing method that when applied 
to a broadband signal, such as that generated by turbulence, acts as a set of band-pass 
filters. This process was applied to data from time-resolved, particle image velocimetry 
measurements of subsonic jets prior to computing the second-order, two-point, space-time 
correlations from which turbulent phase velocities and length and time scales could be 
determined. The application of this method to large sets of simultaneous time histories is 
new. In this initial study, the results are relevant to acoustic analogy source models for 
jet noise prediction. The high frequency portion of the results could provide the turbulent 
values for subgrid scale models for noise that is missed in large-eddy simulations. The 
results are also used to infer that the cross-correlations between different components of the 
decomposed signals at two points in space, neglected in this initial study, are important. 
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Introduction 

a Provide statistical properties of turbulence to aid in 
modeling jet noise sources 
9 Acoustic analogy methods and subgrid models 
® Two-point, space-time correlations 

a numerically - direct numerical simulation 
® experimentally - time-resolved PIV 

a Method - empirical mode decomposition 

9 Automatic filtering for small scales - subgrid models 

a Presentation: 

Q Time-resolved PIV measurements 
Q Empirical mode decomposition - new application 
Q Initial analysis of and results from correlation calculations 
Q Comments on method and correlations 
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Overview of Time-Resolved PIV 


High sampling rate means small field of view 

25 kHz 

~ 1 sec, 24993 points 
178.85 mm by 5.18 mm 
70 by 5 points 
Ax/D = 0.0510 
Ar/D = 0.0255 

Flow conditions for converging nozzle D = 50.8 mm 
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Mean Flow 

SP7 jet with Mj = 0.98 


Axial velocity 
contours 

Axial velocity 


Turbulence 

intensities 



x/D 


Reference points: r/D = -0.49 

x/D = 3.41, 7.82, 10.0 
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Empirical Mode Decomposition (EMD) 


• Adaptive signal processing method for a general, 
non-stationary and nonlinear signal 

a Huang, et al.: Proc. Roy. Society London, 1998 
a Separates signal into basis functions (a posteriori, based 
on the signal, not a priori, based on harmonics) 
a Defined by an algorithm, limited mathematical foundation 
a Applied in vibrations, geophysics, biology, cosmology, 
finance, others 

a Application to large set of simultaneous time histories all 
having the same broadband nature and then computing 
correlations is new 
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EMD Algorithm I 


Start with the signal: the axial velocity fluctuation after 
subtracting the mean axial velocity 



u(i) 
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EMD Algorithm I 


Identify all maxima of the signal 



u(i) 
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EMD Algorithm I 


Interpolating function found using cubic spline 
gives maximum of envelope 
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EMD Algorithm I 


Repeat for the minimum envelope 
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EMD Algorithm I 


Compute average of the envelope 



u(i) 


(i), e . (i) 

x' min' ' 


(e . + e )/2 

' min max' 
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EMD Algorithm II 


An interpretation of EMD is that it separates the local (high 
frequency) oscillations from the local (low frequency) trend 



Qf 
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EMD Algorithm II 


An interpretation of EMD is that it separates the local (high 
frequency) oscillations from the local (low frequency) trend 

signal local trend 



Qf 
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EMD Algorithm II 


Subtract local trend from signal 
Ideally should get zero mean oscillation 

Sift 1 local trend 
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EMD Algorithm II 


Subtract local trend from signal 
Keep sifting 


Sift 2 local trend 
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EMD Algorithm II 


Subtract local trend from signal 
Keep sifting 


Sift 3 local trend 
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EMD Algorithm II 


Subtract local trend from signal 
Keep sifting 


Sift 4 local trend 
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EMD Algorithm II 


Subtract local trend from signal 
Keep sifting 


Sift 5 


local trend 
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EMD Algorithm II 


Sift 6 



Subtract local trend from signal 
Keep sifting 


local trend 
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EMD Algorithm II 


Subtract local trend from signal 
Keep sifting 


Sift 7 


local trend 
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EMD Algorithm II 


Subtract local trend from signal 
Keep sifting 


Sift 8 local trend 



18th AIAA/CEAS AEROACOUSTICS CONFERENCE 


EMD Algorithm II 


Subtract local trend from signal 
Keep sifting 


Sift 9 


local trend 
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EMD Algorithm II 


Subtract local trend from signal 
Stop sifting 


Sift 10 local trend 



Intrinsic Mode Function (IMF) Process with EMD 

- |(extrema) - (zeroes)| < 1 to get next IMF 

- Mean about zero 

- Symmetric 
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Example Intrinsic Mode Functions (IMF) 

0.06 second time sample of a 1 second TR-PIV time history 


Signal 

IMF 1 


IMF 2 


IMF 4 
IMF 6 
IMF 8 
IMF 10 


IMF 12 


4H**fr**1y^ 




. — A/y\^A^AA/w^A/yyw^wy^ — ^\z\^/NY^/^^\/^^yyv'A/y/wv! 



IMF 3 
IMF 5 
IMF 7 
IMF 9 


IMF 11 


Res. 


0.42 


0.44 0.46 

t (sec.) 


0.48 


IMFs nearly orthogonal and uncorrelated: u(t) = Yh n = i C„(f) + z a/(0 
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Broadband Signal Processed by EMD 

EMD acts as a dyadic filter bank 
« Set of overlapping band-pass filters 
« half or double range of neighbor 
o approx, constant shape on log scale 
9 mean frequencies ffj « f 0 2~ n for n-th IMF 



C 

CL) 



Dashed line: 

Add IMFs 2 and 3 


fD/U 


Ref. Point: x/D = 3.41 shear layer 
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Broadband Signal Processed by EMD 

EMD acts as a dyadic filter bank 
« Set of overlapping band-pass filters 
« half or double range of neighbor 
o approx, constant shape on log scale 
9 mean frequencies ffj « f 0 2~ n for n-th IMF 



Dashed line: 

Add IMFs 2 and 3 


fD/U 


Ref. Point: x/D = 7.82 end of potential core 
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Broadband Signal Processed by EMD 


EMD acts as a dyadic filter bank 
« Set of overlapping band-pass filters 
« half or double range of neighbor 
o approx, constant shape on log scale 
9 mean frequencies ffj « f 0 2~ n for n-th IMF 
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tUj/D 


Correlation I 


Total r u {x,w,r) 


u\ (x, Q • u\ {x + m , t + r) 
/ {X, t ) 2 • u\ {x + m f + t ) 2 




Correlation I 



Correlation I 




tUj/D 


Correlation I 




tU/D 


Correlation I 


IMF 4 


rnn{x,rh,T) 


C-\ n{ x i t) ■ C-\ n (x + r]-\,t + r) 

sj U\ ( X , t ) 2 l l\(x + ri^t+^f 



scale: -0.1 to 1 


Ref. Pt.: 
x/D = 7.82 
r/D = -0.49 

Issue within 1 
grid spacing 


'l/D 
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Integral Length and Time Scales 


+oo 

rii(*,»7i ,r = 0)d 77! 


h i (x, m = U c t, t) dr 


-2 - 1.5 -1 - 0.5 0 0.5 1 1.5 

11,/D 

Estimate the integral length and time scales by determining how 
far it takes the correlation to decay by 1/e from the peak value. 
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Integral Length and Time Scales 



Estimate the integral length and time scales by determining how 
far it takes the correlation to decay by 1/e from the peak value. 
Contour resembles an ellipse. 
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& 

W 11 


t u/d 


Ellipse 

Fit an ellipse to the contour where the level is 1/e times the 
peak value of the correlation. Find coefficients a to f. 


(m N 

\ 2 | b ^ tU ^ | cl 

/ tUj\ 

VD/ 

> 1 °D D ' 1 

, D ) 




jUj 

'~D~ 


f = 0 
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Correlation I 



.4 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

-2 -1.5 -1 -0.5 0 0.5 1 1.5 


11/D 


Estimate: 

r+00 

L n = / r u (x,r)i,T = 0 ) d 77 ! 
Jo 

r+00 

r , 7 = / /'ll (X, 771 = U c t , t ) dr 
Jo 

using ellipse eq. 


/-,/D 

r ri Uj/D 


Total 

0.596 

0.366 

3.057 
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Correlation II 



Estimate: 

/»+ OO 


L'n — 


/•“TOO 

/ /'ll m,T = 0) d77i 

Jo 


r+oo 


r 1 1 (x , 771 


U c t, t) dr 


using ellipse eq. 


U c /(/j 

L v /D 

r v Uj/D 


Total IMF 1 
0.596 0.670 

0.366 0.076 

3.057 0.403 
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Correlation II 



-2 - 1.5 -1 - 0.5 0 0.5 

'l/D 

1 1.5 



Total 

IMF 1 

IMF 2 

Uc/Uj 

0.596 

0.670 

0.586 

LJD 

0.366 

0.076 

0.206 

t v Uj/D 

3.057 

0.403 

1.036 


= 0 ) drn 


l c r, t) dr 
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0 


Correlation II 


Q 

3 

H 



Estimate: 


L v 



= 0)d 7?! 



r 1 1 { x , 771 


U c t, t) dr 


using ellipse eq. 


Uc/Uj 

LJD 

r v Uj/D 


Total 

IMF 1 

IMF 2 

IMF 3 

0.596 

0.670 

0.586 

0.540 

0.366 

0.076 

0.206 

0.375 

3.057 

0.403 

1.036 

2.246 




18th AIAA/CEAS AEROACOUSTICS CONFERENCE 


Correlation II 

Estimate: 

r+oo 

Lr, = / r U (X,T)l,T= 0) d77! 
Jo 

r+oo 

r, 7 = / rn (Jf, 771 = L/ C r,r)dr 
Jo 

using ellipse eq. 


-2 - 1.5 -1 - 0.5 0 0.5 1 1.5 

11/D 



Total 

IMF 1 

IMF 2 

IMF 3 

IMF 4 

Uc/Uj 

0.596 

0.670 

0.586 

0.540 

0.384 

LJD 

0.366 

0.076 

0.206 

0.375 

0.545 

t v Uj/D 

3.057 

0.403 

1.036 

2.246 

2.266 



High frequency, small-scale for subgrid model 
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Frequency Dependent Values 


Compute complex coherence 7 ,y (Kerherve et al. ( 2006 )) 


Timescale t ^( x , u ) = / |7 11 (x,? ?1 ,u;)| 

^cl (w) Jo 

Phase velocity u c1 (w) = uj/\d4>(r]-\ , oj)/drj-\ \ 

a Scales difficult to compute without modeling, ex. 3 ? {711} 
oscillates, integration difficult 

a Estimate length scale from (7^ | decays by 1 /e 
a Estimate time scale from I711 |/t/ c1 decays by 1 /e 


Length scale 
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Frequency Dependent Phase Velocity 

Issues: Higher IMFs (lower frequency) lack resolution, low 
number of averages 



Frequency Dependent Phase Velocity 

Compared to ellipse eq. based phase velocity 



Frequency Dependent Phase Velocity 

Compared to ellipse eq. based phase velocity 
Compared to Morris & Zaman (201 0) fit, Mj = 0.25 jet 



Frequency Dependent Phase Velocity 

Compared to ellipse eq. based phase velocity 
Compared to Morris & Zaman (201 0) fit, Mj = 0.25 jet 



Ref. Pt.: 
x/D = 3.41 
r/D = -0.49 
shear layer 
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Frequency Dependent Phase Velocity 

Compared to ellipse eq. based phase velocity 
Compared to Morris & Zaman (201 0) fit, Mj = 0.25 jet 
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Ref. Pt.: 
x/D = 10.0 
r/D = -0.49 
downstream 
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Frequency Dependent Length Scale 


« Compared to ellipse eq. based length scale 

« Using |7n| instead of 3 ? {711} over-estimates the length 
scale 
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Ref. Pt.: 
x/D = 10.0 
r/D = -0.49 
downstream 
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Frequency Dependent Time Scale 

« Compared to ellipse eq. based time scale (solid) 

« Direct interpolated value from correlation function (dashed) 
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Ref. Pt.: 
x/D = 10.0 
r/D = -0.49 
downstream 
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IMF Orthogonality and Correlation 

Accepted that IMFs nearly orthogonal and uncorrelated at a point 


r-ii (*, Vi, T ) 


U\ (x, t) ■ u\ (x + m , t + r) 

\J u\ (X, 0 2 ■ u; (X + ^ , t + t ) 2 


Cln(x. f) • C| n (x + Til , t + t) 
n=1 U'(X, t) z - u\(x + r/i,t + T) z 
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IMF Orthogonality and Correlation 


What about the correlation between two points? 

(x, f) ■ uj (x + 771 , f + r) Ci n (x, t) ■ C| n (x + 771 , f + r) 

n 1 (X, rj-i , rj = = > 

V U 1 (*> f ) 2 ■ u; (X + ^ , f + t) 2 n=1 y U' (X, f) 2 ■ u; (X + ) ?1 , f + r) 2 

+ cross-correlation terms 
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Concluding Remarks 


a New application of empirical mode decomposition 

a equivalent to passing array of time histories through a bank 
of band-pass filters - broadband signals only 
a highest frequency range turbulent statistics may be 
applicable to subgrid modeling 
a Issues: 

a higher sampling rates: higher Strouhal numbers, finer time 
resolution at high frequencies 
a longer time histories: greater frequency resolution, 
increased number of averages 
a finer spatial resolution: capture high frequency changes in 
a short distance 
a Further analysis: 

a IMF cross-correlation terms 
a Other approaches to using EMD results 
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